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Abstract. We present Ha and [Sn] observations for the 
Local Group irregular galaxy IC 1613 using the PUMA 
scanning Fabry-Perot interferometer. Our goal is to an- 
alyze the kinematics of the ionized gas in the complex 
sample of superbubbles located in the whole extension of 
our field (10'), which includes most of the optical emission 
of this galaxy, and to study the inter-relationship between 
young stellar associations and nebulae based on a previ- 
ous study that we have made on the stellar associations of 
the central region of this galaxy. The ionized gas in this 
galaxy is distributed in classical H II regions and in a series 
of superbubbles (also called giant shells) covering a large 
fraction of the optical extent of the galaxy. 

We present a catalog of kinematical properties of both 
the Hn regions of this galaxy and the superbubbles. We 
have also compared the kinematics of the ionized gas in 
H II regions to search for possible dynamic differences be- 
tween neutral and ionized gas. 

Key words: Galaxies: irregular - Galaxies: individual: 
IC 1613 - Galaxies: kinematics and dynamics - Local 
Group - (ISM): HII regions - ISM: bubbles 



1. Introduction 

Irregular galaxies constitute a substantial fraction of all 
galaxies. They are also dominant among the actively star- 
forming galaxies. Because they lack spiral density waves, 
they are excellent laboratories to examine the star for- 
mation process in alternative scenarios. Also, they can 
provide important clues about the mechanisms of forma- 
tion of "bubbles" (the ring-shaped nebulae with diame- 
ters, D ^ 50 pc), "superbubbles" or "giant shells" (ring- 
shaped structures with D 400 pc) and "supershells" 
(ring-shaped structures with D « 1 kpc). Massive stars 
can dramatically change the interstellar medium (ISM) 
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of these galaxies by means of supernova explosions and 
stellar winds. On the other hand, the ISM content and 
energetics can have important effects on the subsequent 
mechanisms of star formation. Indeed, the compressions 
produced by shocks formed by supernova explosions or 
powerful stellar winds on the molecular clouds could trig- 
ger the formation of a second generation of stars at the 
edges of the shells formed by these events (self-induced 
star formation). 

The Local Group irregular galaxies are exceptional tar- 
gets for the study of the inter-relationship between gas 
and stars because they are easily resolved. This means 
that we could know the stellar content of nebular Hn 
complexes and superbubbles and try to identify for rela- 
tionships between the interior stellar associations and the 
dimensions and kinematics of the superbubbles. This kind 
of investigation has been performed for the most favor- 
able case, the Large Magellanic Cloud (Rosado 1986; Oey 
1996a, 1996b, amongst other references). Indeed, it is in 
this galaxy that a clear example of self-induced star for- 
mation has been found: the nebular complex Nil (Parker 
et al. 1996; Rosado et al. 1996). 

In this work, we study the nebulae and massive stel- 
lar content of the nearby irregular galaxy: IC 1613. This 
galaxy belongs to the Local Group and it is located at 
a distance of 725 kpc (Freedman 1988a, 1988b). Table 1 
reviews some of its main properties. Because of its prox- 
imity, it has been the target of several studies. Its stel- 
lar content has been analyzed by Baade (1963), Sandage 
(1971), Sandage & Katem (1976), Hodge (1978, 1980), 
Freedman (1988a, 1988b) and Georgiev et al. (1999; here- 
after Paper I). Several candidate Wolf-Rayet (WR) stars 
were detected but only one of them was spectroscopically 
confirmed over the whole extent of this galaxy (Azzopardi 
et al. 1988; d'Odorico & Rosa 1982; Kingsburgh & Bar- 
low 1995). This WR is a WO-type star embedded in a 
large nebular complex located at the southern region of 
this galaxy (Rosado et al. 2000). The Hn regions of IC 
1613 have been studied by Sandage (1971), Hodge et al. 
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(1990; hereafter HLG), Price et al. (1990) and Hunter et 
al. (1993). The peculiarity of the kinematics of some of the 
superbubbles of ionized gas has been analyzed by Meaburn 
et al. (1988). Their work focused only on the NE region 
of the galaxy where bright H n complexes are easily de- 
tected. Other spectroscopic studies on the nebulae of IC 
1613 mainly examine the only SNR (d'Odorico et al. 1980; 
Pcimbert et al. 1988; Dickel et al. 1985; Goss & Lozinskaya 
1995; Rosado et al. 2000). The distribution of the Hi gas 
in IC 1613 has been analyzed by Lake & Skillman (1989). 
These authors have found that the Hi content makes up 
to 20 percent of the total mass (see also Huchtmeier et al. 
1981). 

In this work, the results of our kinematic analy- 
sis, making use of a scanning Fabry-Perot interferometer 
(PUMA, Rosado et al. 1995), are presented. We analyze 
the relationship between nebulae and the stellar associa- 
tions studied in Paper I. The observations and data reduc- 
tions are described in Section 2. In the remaining sections 
we study the kinematic behavior of the H II regions and 
superbubbles. In Section 3 we show the nebulosities of this 
galaxy in several emission lines. We also indicate how our 
photometric calibration has been done. In Section 4 we 
study the kinematics of HII regions and superbubbles in 
IC 1613. We obtain the peak velocities and velocity disper- 
sions of the H II regions, as well as their stellar content. We 
also obtain the H II luminosity function and diameter dis- 
tribution. In Section 5 we study the superbubbles found, 
obtaining their dimensions, expansion velocities, Ha fluxes 
and studying their stellar content. We also obtain some 
important physical parameters of the superbubbles and 
the superbubble diameter distributions. We analyse the 
possibility that the superbubbles were formed by the com- 
bined action of stellar winds and supernovae explosions. 
In Section 6 we present the radial velocity field in our field 
of view, obtained from diffuse ionized gas, H n regions and 
superbubbles. We compare our trends with the results of 
Hi observations. In Section 7 we present our conclusions. 

2. Observations and data reduction 

The observations were carried out during the nights of 
December 5 and 6, 1996 and November 15, 1998 with the 
UNAM Scanning Fabry-Perot interferometer PUMA at- 
tached to the f/7.9 Ritchey-Chretien focus of the 2.1 m 
telescope at the Observatorio Astronomico Nacional at 
San Pedro Martir, B.C., Mexico. 

The PUMA setup was composed of a scanning Fabry- 
Perot interferometer (SFPI), a focal reducer with a f/3.95 
camera, a filter wheel, a calibration system and a Tek- 
tronix lKxlK CCD detector (Rosado et al. 1995). Table 2 
reports the main characteristics of the instrumental setup. 

We obtained a series of direct images in the lines of Ha, 
[Sil] (A 6717 and 6731 A) and [Om] (A 5007 A) using the 
PUMA in its direct image mode (without the SFPI). The 
exposure times of each of the direct images were 60 s, 60 



s and 300 s for the Ha, [S n] and [O in] , respectively . The 
CCD reading was binned 2x2, giving a pixel size of 1'.'17 
(equivalent to 4.21 pc at the adopted distance) and a field 
of view covering 10' including the NE region and several 
other interesting regions of IC 1613. 

In addition, we obtained scanning FP interferomet- 
ric data cubes in the lines of Ha and [Sn] (A 6717 and 
6731 A). The Ha cubes share the same scale and field of 
view as the direct images. Given that the [S II] intensities 
were lower than the Ha intensities, we binned 4x4 the 
CCD reading of the [Sil] cubes and correspondingly, we 
obtained a pixel size of 2'.'34. 

The phase calibration was made by taking data cubes 
of a calibration lamp before and after the galaxy exposure. 
We used a hydrogen lamp for calibrating the Ha nebular 
cubes and a neon lamp for calibrating the [S n] nebular 
cubes. Since the redshift of this galaxy is not high, phase 
corrections, due to the difference in wavelength between 
the galaxy and the calibration line, are not required. 

The FP data cubes are composed of 48 steps, with inte- 
gration times of 60 seconds per step. In this form, the final 
data cube dimensions in Ha and [Sil] were 512x512x48 
and 256x256x48, respectively. For the Ha observations 
we obtained 4 data cubes that were summed giving a to- 
tal exposure time of 3.2 hours. For the [S n] observations 
we obtained 3 data cubes that also were summed, with a 
total exposure time of 2.4 hours. 

Reductions were carried out using the software pack- 
age CIGALE (Le Coarer et al. 1993). Our data cubes in 
Ha are contaminated with line-sky emission such as the 
geocoronal Ha line and the lines of OH at 6577.18 A and 
6577.38 A, respectively. The night-sky lines fell well out- 
side the main emission of IC 1613, due to the systemic ve- 
locity of this galaxy. These night-sky lines were subtracted 
using an interactive routine of CIGALE. The data cubes 
in the [S II] lines are less contaminated by line-sky emis- 
sion; however, the S/N ratio is not as good as that of the 
Ha cubes and so will be used only for comparison. 

To extract the kinematic information from the FP data 
cubes we proceeded in different ways for the H II regions 
and for the superbubbles. In the case of the H n regions, we 
obtained only one radial velocity profile at Ha, integrated 
over a box containing the whole Hn region. In the case of 
the superbubbles, we proceed in a different way in order to 
seek for internal motions revealed by line splitting of the 
velocity profiles. We obtained the radial velocity field of 
the superbubble by dividing the superbubble dimensions 
into several small boxes over which we obtained several 
radial velocity profiles, both in Ha and [S II] . 

The radial velocity profiles were fitted by Gaussian 
functions after deconvolution by the instrumental func- 
tion (an Airy function). We constructed a map of the 
monochromatic (continuum) emission by integration of 
the radial velocity profile of each pixel, up (down) to a 
certain fraction of the peak. This map enabled us to sepa- 
rate the monochromatic emission from the continuum (Le 
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Coarer et al. 1993). To get the radial velocity map of our 
observed field, we found the central velocity in the radial 
velocity profile of each pixel, that has a level above 2 a 
the value of the standard background emission, to ensure 
galaxy membership. 

3. Morphology of the nebulae in IC 1613 

Figure 1 shows the emission of IC 1613 in the lines of 
Ha, [S n] and [Om] as obtained from our direct im- 
ages. The photometrical center of IC 1613 is located at 
RA(1950)=01:02:19.6 and Dec(1950)=+01:51:56.1 (Ma- 
teo 1998). The most active region of IC 1613 is located 3' 
(or 663 pc) North-East from the photometric center. Here- 
after, we will call this quadrant the NE region. Bright Hil 
regions, bubbles and supcrbubbles arc situated in the NE 
region. Numerous H n regions and larger structures, some 
of them ring-shaped, are superimposed over the extent 
of the chaotic body of the galaxy. Except for the bright 
network of supcrbubbles located in the NE region of the 
galaxy, the nebulae are weaker in the [S il] lines and they 
do not show any appreciable [O in] emission. 

In the South of this rich sector of the galaxy (at 3f8 
or 838 pc from the center), there is a giant Hn region 
(Sandagc 3). Southwest of the center of the galaxy at 2f4 
(535 pc) we located the Hn region Sandagc 2. 

3.1. The Hn region content 

The Hn regions of IC 1613 have been catalogued by 
Sandage (1971), Lequeux et al. (1987), Price et al. (1990) 
and HLG (who catalogued 77 H n regions, reporting their 
diameters and Ha fluxes). 

Our observed field contains 44 diffuse nebulae accord- 
ing to HLG. The limits of each Hn region or group of 
H II regions were defined by the contour where the inten- 
sity of the Ha emission falls to the average local intensity 
of the diffuse background. However, in some cases, the 
boundaries between diffuse Hn regions are not obvious. 
In addition, we were not able to resolve some regions. In 
these cases, we assign to them a common boundary, and, 
consequently, only 30 regions were clearly defined. 

3.2. Photometric calibration 

We calibrate our monochromatic Ha image in flux, using 
the measurements of HLG. Eight isolated and round re- 
gions were selected from HLG lists (Table 3). For each 
region we measured the flux of our monochromatic Ha 
image, using apertures with the same equivalent area as 
in HLG. A least squares fit to the data gives a conver- 
sion factor of 1 count s _1 (in our Ha image) = 6.26 10~ 18 
erg cm~ 2 s _1 . We preferred a conversion factor per unit 
area for easy calculations of the flux in irregular HII re- 
gions and supcrbubbles. Polygonal apertures were chosen 
such that they must contain all emission from a given H II 



Fig. 1. Top to bottom: Direct images of IC 1613 in Ha, [S n] 
and [Om], respectively. The exposure time is 60 sec for Ha 
and [Sn], while for [O ill] it is 300 sec (see text for details). 
The scale and the orientation are indicated in the upper panel. 

region or superbubble. Pixels within the apertures were 
summed to generate the total Hn region or superbub- 
ble plus background counts. The background emission was 
computed using apertures with the same area as the ob- 
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ject, as close to the object as possible. The total Ha flux 
for every H n region and every superbubble was obtained 
when the background was subtracted and the resulting 
counts multiplied by the conversion factor. Ha fluxes are 
presented in column 9 of Table 4 for the H n regions and 
column 6 of Table 6 for the superbubbles. Figure 2 shows 
the comparison between our measurements and HLG for 
the common HII regions not included as calibrators. 




10 13 2x10 13 3x10 13 

This work Flux [erg cm 2 s 1 ] 



Fig. 2. Comparison between our flux determinations and HLG 
fluxes. The calibrator HII regions were not used in the fitting. 
The values of the correlation coefficient (r) and the slope (b) 
are given to the left. The dotted line is our best fit, and the 
solid line is the diagonal. 



4. Kinematics of the Hn regions in IC 1613 

Figure 3 shows the Ha monochromatic image of IC 1613 
(continuum subtracted) obtained from our FP cubes us- 
ing the method described in Section 2. Superimposed on 
this image, we marked the rectangular boxes over which 
the radial velocity measurements of the observed Hn re- 
gions were performed. We obtained radial velocity pro- 
files of these sources by integrating in those rectangular 
boxes. We will refer to these profiles as "integrated" be- 
cause they are obtained by integrating over the whole ex- 
tent of the Hn region. Also, in this figure, we marked with 
asterisks the WR star candidates (Armandroff & Massey 
1985). Only one of these candidates was confirmed as a 
WR star; another was confirmed as the SNR Sandage 8 
and the remaining were identified as OB stars (Azzopardi 
et al. 1988). 

The integrated radial velocity profiles are, in general, 
single for all the H n regions; in this case, a single Gaussian 
function was fitted as described above. In only two cases: 
the SNR Sandage 8 (our source N18, see Table 4) and 
the giant shell Sandage 10/13 (N27), the integrated radial 
velocity profiles were complex. In those cases, we carried 



out a profile decomposition into two Gaussian functions. 
We then list the two velocity components obtained from 
our profile fitting. A more detailed study of the kinematics 
of S8 from the FP data will be given elsewhere (Rosado 
et al. 2000). Also, a more detailed study of the internal 
kinematics of the superbubble N27 = 67 (HLG) and of 
other superbubbles in this galaxy is given in the following 
section. 

Table 4 gives a catalogue of the kinematical parame- 
ters of the H II regions of IC 1613, within the field of view 
of our observations. It is important to note that some of 
the objects listed in this table are ring-shaped and conse- 
quently are better classified as superbubbles and will be 
studied in more detail in the following section. We include 
these objects in Table 4 because they appear in previous 
catalogues (Sandage 1971 and HLG) but we exclude them 
from the calculation of the H n region luminosity function 
and diameter distributions. 

In Table 4, we have listed the name (adopting both 
HLG's and Sandage's conventions), diameter (see the def- 
inition for this parameter in Section 4.2), peak velocity 
and velocity dispersion from the Gaussian fitting to the 
radial velocity profile, and the Ha flux, derived as de- 
scribed in Section 3.2. To obtain luminosities, our fluxes 
were corrected only for galactic foreground extinction us- 
ing the Ab=0.21 mag value given in the RC2 catalog (de 
Vaucouleurs et al. 1976). The effects of internal extinction 
in these determinations seems to be not crucial (Hippelein 
1986). The radial velocities were corrected to the heliocen- 
tric system of reference. Our dispersions {(Jobs) have been 
corrected for instrumental (<7i ns t), thermal {(Jth) and in- 
trinsic {(Jintrin) broadenings respectively, according to: 

2 _ 2 _ 2 _ 2 _ 2 

corr obs °inst °th intrin 

where <Ji nst — 43.5 km s _1 , a t h = 9.1 km s _1 for hydrogen 
gas at T e =10 4 (Spitzer 1978) and (Jintrin — 3 km s _1 for 
the fine structure of the Ha line. Some authors also refer 
to cr corr as the turbulent velocity dispersion. 

An inspection of this table shows that the heliocen- 
tric radial velocities range from -228 to -254 km s _1 . The 
equivalent diameters range from 27 pc (corresponding to 
region N24) to 253.6 pc (for region Sandage 3 or N9) and 
their mean value is around 82 pc. The corrected velocity 
dispersions are supersonic for most of the objects with val- 
ues ranging from 16 to 35 km s -1 with a mean value of 20 
km s _1 . Some Hn regions have low S/N ratios and conse- 
quently their determinations of peak velocity and velocity 
dispersion are uncertain. These cases have been quoted 
within parentheses in Table 4 . The Ha fluxes range from 
1.60xl0~ 15 to 9.24xl0~ 13 ergs cm~ 2 s _1 ; these lower and 
upper limits correspond to the H n regions N13 and N27, 
respectively. 

In Table 5 we give the stellar content in order to iden- 
tify the exciting star or stars that are related to the ana- 
lyzed H II regions. Column 1 gives the identification labels 
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Fig. 3. Ha monochromatic image of IC 1613. The background has been subtracted at 2a its mean value in order to reveal the 
weaker features. The boxes correspond to places where the kinematical analysis was performed. Sandage identifications have 
been quoted between parentheses if relevant. The asterisks indicate positions for the eight WR star candidates (Armandroff & 
Massey 1985). Some labels are denoted with crosses, at the left bottom of the figure. The scale and orientation are indicated to 
the left and right, respectively. In this Figure the field covers approximately 9' x 9'. 



used in this work. Column 2 and 3 list the WR candidates 
(Armandroff & Massey 1985) that coincide with some of 
our H II regions and their spectral type (or related object) 
according to Azzopardi et al. (1988). Columns 4, 6 and 7 
indicate the stellar associations and clusters determined in 
previous works. All but four H n regions have an evident 
stellar content that is probably the source of the UV flux. 

4-1. The Hil region luminosity function 

Figure 4 shows the Ha luminosity function for the H n re- 
gions in IC 1613. For several galaxies, it has been shown 
that the brighter luminosity end of the H II region luminos- 
ity function can be fitted with a power law in the form: 
A^(L)=AL° dL , where N is the number of Hn regions 
with luminosity L in an interval of luminosity dL, and the 



slope a is related to the morphological type of the galaxy 
(Kcnnicutt ct al. 1989), ranging from mean values of -2.3 
(Sab-Sb galaxies) to -1.75 (Irregulars). It also has been 
suggested that different ranges of the luminosity function 
are related to different kind of objects, such as Hil re- 
gions ionized by single stars and those ionized by clusters 
of stars (Hodge et al. 1989; Youngblood & Hunter 1999). 

To obtain the H n region luminosity function we take 
the logarithm of the number of H n regions divided by 
the corresponding luminosity bin width (0.26 dex in our 
calculation). The bin width was chosen in terms of our 
relatively small number of objects (30 H n regions) and the 
fit was calculated to the bins brighter than logL = 36.5. 
Our best fit slope for IC 1613 gives a = -1.54 ± 0.1. This 
value is similar to the value (a = — 1.6 ± 0.2) found by 
HLG, and is within the mean values found for a sample 
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of irregular galaxies studied by Kennicutt ct al. (1989), 
Strobel et al. (1991), Kingsburgh & McCall (1998) and 
Youngblood & Hunter (1999). 



4-2. The Hll region diameter distribution 



We measured the diameter of the Hll regions on the 
monochromatic Ha image and they are reported in col- 
umn 4 in Table 4. The reported diameter is the equiva- 
lent diameter, which is defined as D = (AA/tt) 1 ^ 2 . In this 
definition, A is the area of the apertures wherein the pho- 
tometry was performed. The measured diameters range 
from 27 to 254 pc (this later value corresponds to large 
H II complexes that are not ring-shaped) . 



The cumulative Hn region diameter distribution 
iV(> D) has been calculated and in Figure 5 we plot 
log iV(> D) vs D. This functional dependence was pro- 
posed by van den Bergh (1981) as a universal law for all 
galaxy types and it has been used to study the statistics of 
the Hn region populations in combination with the Hn 
region luminosity function. The slope in the cumulative 
H II region diameter distribution is related to the charac- 
teristic diameter Dq of the Hn regions in a galaxy. For 
IC 1613, a least squares fit gives Dq = 65.41 ± 3 pc. This 
value is marginally in agreement with a previous one de- 
termined by HLG who obtained Dq=56 pc. However our 
value differs from that of Price et al. (1990) who quoted 
a larger value (Do=83 pc). This discrepancy can be easily 
explained, because Price et al. (1990) also considered some 
superbubbles in their calculation (see Section 5). Looking 
for other correlations, Hodge (1983), using a sample of 18 
irregular and related galaxies (IC 1613 was not in his sam- 
ple), found a good correlation between the characteristic 
diameter Dq and the absolute magnitude Mb of a galaxy. 
In a subsequent work, Strobel et al. (1991) increased 
the sample of Hodge (1983) with 7 additional galaxies 
(IC 1613 was also not included in these new objects). 
They found that the previous correlation between Dq 
and the absolute magnitude Mb becomes weaker, having 
the form logA, = 0.25(±0.29) - 0.098(±0.017)M B . For 
IC 1613 this relation gives a value of log£>o = 1.64 - 
taking Mb = —14.2 from Mateo (1998) - in agreement 
with our calculation (log_D = 1-83). However, if we plot 
our resulting Do in Figure 11 in Strobel et al. 1991, our 
data point falls outside the line fit in that plot, a fact 
that points towards the conclusion given by Strobel et al. 
(1991). In that sense, in a more recent work, Youngblood 
& Hunter (1999), using a sample of 29 irregular and Blue 
Compact Dwarfs (IC 1613 not included), have shown that 
this correlation is indeed a very weak one (see their Figure 
9). 




log L 

Fig. 4. The Ha luminosity function for the HII regions in IC 
1613. The dashed line is the least squares fit to a power law. 
The slope value is given to the right. 




D (pc) 

Fig. 5. The cumulative size distribution for the Hll regions in 
IC 1613, where the dashed straight line is a least squares fit. 
The value for the characteristic diameter is given to the right. 



5. Kinematics of the superbubbles (giant shells) 
of IC 1613 

Figure 6 shows one of the velocity maps (at V/ le ;io=-234 
km s _1 ) of the Ha emission of IC 1613, obtained from 
our FP observations. The contrast of this image is higher 
than the direct image shown in Figure 1 because the FP 
interferometer acts as a filter with a very narrow width 
(0.41 A) while the bandwidth of the filter used to obtain 
the direct image in Figure 1 is 20 A. Thus, it allows the 
detection of the line emission from the galaxy while reduc- 
ing spurious line or continuum emission. Furthermore, in 
order to enhance the detection of low emissivity features 
of large dimensions, this velocity map was smoothed spa- 
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tially using a Gaussian filter of 4 pixels width in both X 
and Y axis. 

We note that the entire observed field is covered by 
networks of shells, rings and filaments whose diameters 
range from 110 pc (30"6) to 325 pc (90"6). We will follow 
the terminology of this type of structures given in the 
introduction. In accord with this terminology, Figure 6 
shows that the entire field of our observations of IC 1613 is 
covered by superbubbles. We were able to differentiate 17 
large diameter ring-shaped structures. The rough location 
of these superbubbles are marked on Figure 6. We detect 
the seven superbubbles or giant shells, GS, (from GSla 
and b to GS6) already reported by Meaburn et al. (1988), 
as well as several others located outside the field (to the S 
and to the W) covered by Meaburn et al. (1988). Most of 
the newly detected superbubbles lie outside the bright NE 
region of the galaxy and they are very faint and difficult 
to visualize. Some of them are coincident with the Hll 
regions and filaments reported by Price et al. (1990) and 
Hunter et al. (1993). 

Figure 7 and Figure 8 show some of the Ha and [S n] 
velocity maps (or channels) where the network of super- 
bubbles in the NE region is clearly appreciated. The helio- 
centric radial velocities for these superbubbles range from 
-100 to -350 km s _1 . The superbubbles appear different 
in the several FP velocity maps, allowing us to disentan- 
gle several of them that appear superimposed in the di- 
rect images. In particular, the giant shells la and lb of 
Meaburn et al. (1988) (corresponding to superbubbles Rl 
and R2, respectively) appear in different velocity chan- 
nels: Rl (Meaburn's GS la) is better visualized in chan- 
nels covering heliocentric velocities from -272 to -253 km 
s _1 (and even more negative velocities not shown in the 
figures), whereas R2 (Meaburn's GS lb) is better appre- 
ciated in channels where the heliocentric velocity ranges 
from -234 to -196 kms" 1 . 

In order to quantify the relative importance of the su- 
perbubbles and to compare with other galaxies, we have 
calculated the 2-D porosity parameter, Q 2 d, in our field, 
defined as the ratio of the total area occupied by the cat- 
alogued superbubbles to the total area of the galaxy seen 
in our field of view (Oey & Clarke 1997). We obtain Q 2 d 
= 0.75, a value that is intermediate between the observed 
Q 2D values for the SMC (1.6) and M 31 (0.01) according 
with Oey & Clarke (1997). However, the Q 2 d of the whole 
galaxy must await observations of the type described in 
this work covering the whole dimension of IC 1613. 

According to the procedure discussed at the end of 
Section 2, we obtained the radial velocity field of the su- 
perbubbles. We used mainly the Ha data; however, we 
used also the [S 11] data for the superbubbles which showed 
bright emission in the [S 11] lines. Some of the superbubbles 
of the NE region appear also in Table 4. In that table we 
have only quoted the velocity values extracted from the 
integrated radial velocity profiles, while in the following 



we will use the radial velocity field in seeking for internal 
motions. 

Figure 9 shows, as an example, the radial velocity pro- 
file of one of the superbubbles where the splitting is clearly 
visible. As we can see from this figure, some of the profiles 
are complex and require to be fitted by a superposition of 
at least two Gaussian functions. We interpret the different 
components as being due to an expansion motion. Table 6 
gives a catalogue of the detected superbubbles, their iden- 
tification (if any) according to Meaburn et al. (1988), Price 
et al. (1990) and Hunter et al. (1993, hereafter HHG), their 
angular dimensions, the peak heliocentric velocities of the 
Gaussian functions fitted to the velocity profile and the 
calibrated Ha fluxes and luminosities. The last column of 
Table 6 corresponds to comments on the morphology of 
the superbubbles. We have found that some of the super- 
bubbles have the "ring of Hll regions" structure (i.e., the 
shell has smaller Hll regions at its boundaries). This mor- 
phological type presents some evidence of shock-induced 
star formation as Rosado et al. (1996) have found for su- 
perbubbles in the LMC. 

5.1. The stellar content of the superbubbles 

The superbubbles catalogued in Table 6 are directly re- 
lated to the stellar associations of IC 1613. These stellar 
associations were catalogued by Hodge (1978) for the en- 
tire galaxy. In Paper I, we have studied the stellar con- 
tent of the NE region of the galaxy and applied an auto- 
matic method to obtain the boundaries of the stellar asso- 
ciations. Figure 10 shows the stellar associations studied 
in Paper I, where the superbubble boundaries shown in 
Figure 6 are also marked. We see that there is a good 
agreement between the locations of most of the stellar as- 
sociations and the centers of the superbubbles. For the 
superbubbles outside the NE region we also find a good 
agreement between the positions of Hodge's (1978) asso- 
ciations and the superbubbles. 

The stellar associations within the superbubbles 
boundaries are reported in Table 7. In this table we have 
listed (1) the superbubble's number (according to our cat- 
alogue), (2) the Paper I and Hodge's (1978) stellar asso- 
ciations located within the superbubble boundaries, (3) 
the number of blue stars of the interior associations (stars 
with M v < -2.0, (B - V) < 0.0 and (U - B) < -0.5, 
corresponding to stars with masses larger than 7 Mq , 
(4) the age of the stellar associations, the possible exis- 
tence of (5) OB stars and (6) blue supergiants with known 
spectral types and (7) the T exponent of the Initial Mass 
Function of the interior associations. Columns 3, 4, 6 and 
7 were derived from the data published in Paper I. Col- 
umn 5 is based on the spectra published for several Wolf- 
Rayet (WR) star candidates in this galaxy (Armandroff & 
Massey 1985, Azzopardi et al. 1988). An examination of 
this table shows that the stellar associations reported in 
Paper I have a better correlation with the catalogued su- 
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Fig. 6. Velocity map (at Vh e ji =-234 km s 1 ) of the Ha emission of IC 1613, obtained from the scanning FP observations. 
Superimposed on this map are the boundaries of the 17 superbubbles catalogued in Table 6. 



perbubbles than Hodge's. This is because Hodge's (1978) 
associations are, in general, larger than those reported in 
Paper I. Given that so many of them are interior to the 
superbubbles, we can conclude that the method of identi- 
fication of stellar associations discussed in Paper I works 
very well. We also see that every catalogued superbubble 
can be associated with a relatively young stellar associ- 
ation (except R4, which probably has an older interior 
association), even the dimmer superbubbles, located far 
away from the bright NE region of IC 1613. This suggests 
that the interior associations are responsible for the for- 
mation of the superbubbles. 



5. 2. Physical parameters of the superbubbles 

From the observed data (listed in Table 6) we can extract 
several physical parameters of the superbubbles: 



1. The rms electron density of a superbubble is calcu- 
lated from its Ha flux and dimensions according to the 
relation: 

rj2(rms) (cm" 6 ) = 2.74 x 10 18 F(Ha) RT 1 (T 2 

This relation was obtained by assuming that the emit- 
ting gas is at T = 10 4 K and distributed in an spherical 
shell of radius R (in pc) and thickness, AR = R/12. The 
other quantities, F(Ha) and 9 2 , are listed in Table 6. 

2. The expansion velocity Vexp, of a superbubble can 
be estimated as one half of the difference in velocities of 
the components listed in Table 6. 

3. The mechanical energy of the superbubble, K.E. = 
0.5 M Vf^jj-p, is calculated according with the relation: 

K.E.(ergs) = 1.7xl0 41 n e (rms) (cm~ 3 ) V 2 EXP (km s" 1 ) R 3 

In obtaining this relation we have assumed that the 
mass is concentrated in the spherical shell of radius R (in 
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Fig. 7. Ha velocity maps of the NE region in IC 1613 where the superbubble network appears bright. In this map, the 
superbubbles Rl, R2, R4, R5, R6, R7, R8, R9, RIO and R13 are more clearly visible because of the high contrast of the print. 

pc) and thickness, AR = R/12. We have also assumed 
that \i — 0.65, as for an ionized hydrogen gas. 

4. The kinematic age of the superbubble can be ob- 
tained from the following relation: t(yr) = 4.8 x 10 5 R/Vexp (km s _1 ) 
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Fig. 8. [SII] velocity maps of the NE region in IC 1613 corresponding to the same velocities as in Figure 7. The bright feature 
to the South corresponds to the SNR Sandage 8. 

This relation has been obtained by noting that the Table 8 lists the physical parameters of the catalogued 
kinematic age of wind blown bubbles can be estimated as: superbubbles calculated as discussed above: the supcrbub- 
t = 0.5 - 0.6 R/V (McCray 1977). ble's number (according to our catalogue), the shell radius 

(see Sect. 5.3), the rms electron density, the expansion ve- 
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locity, the mechanical energy and the dynamical age of the 
superbubble. Meaburn et al. (1988) have obtained expan- 
sion velocities for five of these superbubbles. Our results 
agree with these previous observations except for the di- 
mensions of some of the superbubbles. 



5.3. The diameter distribution of the superbubbles 

The angular dimensions of the superbubbles were deter- 
mined on the velocity maps at Ha and they are given in 
column 3 in Table 6. The reported dimensions in this ta- 
ble correspond to the extension of the semiaxes in arc sec. 
In the cases where we have two different values, i.e., when 
the shell is not circular, we average them. The resulting 
linear radii are listed in column 2 in Table 8. 

Comparing the size distribution of the H n regions and 
Hll complexes of IC 1613 (see Section 4.2) against the 
superbubble size distribution we note the following dif- 
ferences: while the diameters of the Hll regions range 
from 27 to 254 pc (this later value corresponds to large 
Hll complexes not ring-shaped), those of the superbub- 
bles range from 110 to 304 pc (note that there is an over- 
lap between these two populations). Thus, the superbub- 
bles have, in general, larger dimensions than those of Hn 
regions, pointing towards a more evolved nebular popu- 
lation. Figure 11 shows the size distribution, N(R) — i?, 
of the superbubbles in IC 1613. It is worth comparing 
this size distribution with the size distributions of other 
galaxies. In order to do so, we have examined the size 
distributions of the galaxies SMC, Holmberg II, M31 and 
M33 reported in Oey & Clarke (1997; see their Figures 
3 and 4). It is important to take into account the small 
number statistics of the superbubbles in the studied field 
of IC 1613 (only 17). Furthermore, our superbubbles are 
ionized regions while the superbubbles in Oey & Clarke 
(1997) correspond to Hi holes. Nevertheless, we see that 
the distributions show some similarity. The peak of maxi- 
mum frequency of the superbubbles in these galaxies cor- 
responds to D « 80 pc (SMC), 120 pc (M33), 180 pc 
(M31) and 500 pc (Hoi II); this peak, in the case of IC 
1613, is at « 100 — 140 pc, thus, similar to the one of M 
33 H I holes. On the other hand, there is a difference in the 
maximum linear extent of the superbubbles: the superbub- 
bles in IC 1613 reach D « 300 pc while the superbubbles 
in the other studied galaxies reach diameters of 720 pc 
(M31), 760 pc (SMC), 1160 pc (M33) and 1400 pc (Hoi 
II); note that the latter two values fall into the category 
of supershclls according to the typical terminology. This 
difference could be due to the fact that we are only detect- 
ing ionized shells while the other studies are based on H I 
holes. We have also obtained the cumulative diameter dis- 
tribution N(> D) for the superbubbles. In Figure 12 we 
plot log7V(> D)vsD. For the superbubbles of IC 1613 
a least squares fit gives Do = 91.84 ± 3 pc. This value 
is in agreement with the value (Dq=83 pc) determined 



by Price et al. (1990) who included some superbubbles in 
their calculation of the diameter distribution. 




Vhel (km/s) 

Fig. 9. Radial velocity profile of one region of superbubble R2 
catalogued in Table 6 and its decomposition in two Gaussian 
functions, (0) and (1) with peak radial velocities of - 244 and 
- 147 km s _1 , respectively, for this particular case. The hori- 
zontal axis plots the heliocentric radial velocity while the ver- 
tical axis corresponds to the intensity (in arbitrary units). The 
low velocity peaks correspond to night-sky lines mixed with 
diffuse emission. 



5.4- Discussion 

We have seen in the preceding sections (more precisely 
in Section 5.1) that every superbubble detected in this 
work corresponds to one (or several) stellar association(s) 
interior to the superbubble. This suggests that the su- 
perbubbles are formed by the combined action of stellar 
winds and supernova explosions of the massive stars of 
the associations. In what follows, we will examine in more 
detail this possibility and the models used to describe this 
formation and evolution. 

The shell formation and evolution as a result of the 
interaction of strong winds of a massive star with its am- 
bient medium has been studied a long time ago by sev- 
eral authors (Dyson & de Vries 1972; Castor et al. 1975; 
Steigman, Strittmatter & Williams 1975; Weaver et al. 
1977, amongst others). The more favored (and complete) 
description of this interaction is the "standard model" 
(Weaver et al. 1977) which considers that a star at rest 
releases a constant wind power, Lw, that interacts with 
the ambient medium (assumed to be homogeneous). As 
a result of this interaction, two shocks are formed (one 
in the stellar wind and another in the ambient medium). 
The shell formation is the result of these shocks and the 
shell expansion is driven by the thermal pressure of the 
shocked stellar wind. The shocked stellar wind has no sig- 
nificant radiative cooling and for this reason this model 
is also called adiabatic or energy conserving. This model 
could also be applied to explain the formation of larger 
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Fig. 10. Sketch of the stellar associations studied in Paper I (shadowed regions). We have superimposed on this some of the 
superbubble boundaries already shown in Figure 6. The scale and the orientation are indicated down to the left and top to the 
right, respectively. 

- The Ha luminosities of the superbubbles vary from 
2.11 x 10 36 to 5.89 x 10 37 erg s _1 . These luminosity val- 
ues imply a number of ionizing photons s" 1 ranging from 
1.55 x 10 48 to 4.33 x 10 49 . These ionizing photons can 
be easily supplied by the interior blue stars. Thus, the su- 
perbubbles seem to be ionized by the stars of the interior 
associations. 

- The dynamical ages of the superbubbles (column 6 
of Table 8) are of the order of 1-2 Myr and in every case, 
shorter than the age of the interior stellar association (col- 
umn 4 of Table 5). The fact that the dynamical ages are 
shorter than those of the stellar associations suggests that 
the stellar associations could form the superbubbles. Fur- 
thermore, the fact that there is a delay (of about 10 Myr) 
between the formation of the stellar association and the 
formation of the superbubble (estimated by its dynami- 
cal age), suggests that only the most massive stars of the 
interior association are responsible for superbubble forma- 
tion when they evolve into supergiants (or perhaps have 
exploded as supernovae). 

- The obtained expansion velocities of the superbub- 
bles are in general quite high. All the known expansion ve- 
locities are larger than 25 km s _1 and there are two super- 
bubbles with expansion velocities larger than 40 km s -1 . 
This is similar to the case of the superbubbles in the LMC 
where Rosado et al. (1981; 1982), Rosado (1986) and Oey 
(1996a) have found several superbubbles with high expan- 



shells (superbubbles and supershells) formed by the com- 
bined action of stellar winds and supernova explosions. 
In spite of the popularity and theoretical simplicity of the 
standard model, several observational difficulties have pre- 
vented its test in real cases. The main difficulty is the poor 
knowledge of the interior stars and of their wind power. 
Observational works (Oey 1996b) where the stellar content 
has been thoughtfully analyzed and using numerical (in- 
stead of the analytical) models that take into account the 
temporal variations of the wind power as a result of the 
stellar evolution (Oey 1996a), have found discrepancies 
between the predicted shell radius and expansion velocity 
in the sense that the expansion velocities are too high for 
the observed radii ( "dynamical discrepancy" ) or between 
the shell radius and velocity and the wind power, L^, 
and pre-shock density, no, in the sense that Lw/ n o seems 
to be overestimated by a factor of 10 ("growth rate dis- 
crepancy"). Having these considerations in mind we will 
proceed to analyze our results, summarized in Table 8. 

A close examination of the physical parameters re- 
ported in Table 8 shows that: 

- The rms electron densities of the superbubbles are 
quite low (from 0.3 to up to 6.3 cm -3 ). This makes it 
quite difficult to obtain electron densities from [S n] line- 
ratio diagnosis because the line-ratios will fall in the low 
density regime (n e < 300 cm ). 
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sion velocities. The kinematics of this class of superbub- 
bles ("high velocity superbubbles" ) cannot be explained 
in terms of the "standard model" of adiabatic wind-blown 
bubbles (Weaver et al. 1977) because the expansion veloc- 
ities are too large for the observed shell radius, according 
to the predictions of this model. 

IC 1613 




R [pc] 

Fig. 11. Size distribution function, N(R) — R, for the super- 
bubbles in IC 1613. 




Fig. 12. The cumulative size distribution for the superbubbles 
in IC 1613. The dashed straight line is a least square fit. The 
value for the characteristic diameter is given to the right. 



In the case of the superbubbles of IC 1613, our stellar 
photometry can give us some knowledge of the massive 
stars interior to the catalogued superbubbles but we can- 



not follow the treatment given in Oey (1996a) because 
further spectroscopic work is needed in order to know ac- 
curately the content of massive stars. In the absense of 
this accurate knowledge we can, nevertheless, estimate the 
required wind power {Lw) of each superbubble and com- 
pare it with the number of massive stars that we detect in 
our photometric observations, in order to see, in terms of 
energy considerations, whether the standard model could 
explain the superbubble formation. 

Assuming that the superbubbles evolve according to 
the standard model of wind-blown bubbles (Weaver et al. 
1977), the wind power, Lw, released by the stars con- 
stantly during all of the dynamical age of the shell is given 
by: 

L w = 3.2xl(r 7 no(cm- 3 )V E xp 3 (km s" 1 ) R 2 xl0 36 ergs s~ : 

where, n , the ambient pre-shock density, is taken as: 
no = n e / 4 by assuming that the swept-up mass in the 
shell comes from a homogeneous sphere of radius R. Both 
n e and R are listed in Table 8. 

A wind power of 4xl0 36 ergs s _1 is reasonably ex- 
pected for one 09.5 or BO supergiant star (Pauldrach et 
al. 1986) during its lifetime in this phase (about 10 6 yr). 
More massive stars can have larger wind powers but the 
duration of those phases is ten times shorter. That is the 
case for 041 stars that can release more than 100 times 
the wind power of lxlO 36 ergs s _1 , or of WR stars that 
can release about 50xl0 36 ergs s _1 ; however, these events 
occur only during time scales shorter than the dynami- 
cal ages of the superbubbles considered here (a few 10 5 yr 
compared with the dynamical ages of about 1 Myr derived 
in this work). 

In Table 9 we list the required wind power Lw that 
can explain the observed expansion velocities and radii of 
the superbubbles studied here, in terms of the standard 
model, and the number of detected blue stars (with My 
< -2, or M* > 7 M ), Nb, quoted in column 3 of Table 7. 
We see that superbubbles with expansion velocities lower 
than 30 km s _1 have wind powers that are easily explained 
by the number of blue stars inside them (assuming that 
they are all supergiants and earlier than BO) . On the con- 
trary, there are some extreme cases, such as R2, where the 
required wind power seems to imply that some of the inte- 
rior blue stars should release much more powerful winds. 
Given our limitations in the accurate knowledge of the 
stellar content of the superbubbles, we can only state that, 
from energy considerations, the superbubbles studied here 
could be formed by the stellar winds of the interior stars. 

It is interesting to note from Table 6 that the su- 
perbubbles of IC 1613 with higher velocities seem to be 
brighter in [S n] emission. On the other hand, an inspec- 
tion of Table 8 shows that all the superbubbles have me- 
chanical energies of at most 80 % the energy of one su- 
pernova explosion (approximately l-3xl0 51 ergs). Taking 
into account that between 20 to 30 % of this energy is 
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transmitted to the ambient medium as kinetic energy, one 
or two SN explosions could also account for the formation 
(or acceleration) of the superbubbles. However, no traces 
of these SN explosions (such as extended X-rays sources 
or non-thermal radio emission) have been found. 

6. The velocity field of IC 1613 

Figure 13 shows the global radial velocity field (not cor- 
rected for projection effects) of the central region of IC 
1613 obtained from our FP data in Ha, considering dif- 
fuse gas, Hn regions and superbubbles. In this map a 
chaotic distribution of velocities across the galaxy can be 
seen: some local events are superimposed onto the regu- 
lar field. These local events are the peculiar kinematical 
signatures of the superbubbles, which are notorious. How- 
ever, when the expansions are extracted, we see a variation 
in the mean velocities across the field. Along the NE re- 
gion, a gradient of velocities ranges from -210 km s _1 in 
the southeast to -250 km s _1 in the northwest. The ex- 
treme detected velocities are due to the contamination of 
the superbubble and SNR internal motions. Indeed, the 
region N27 (with the maximum velocity of —139 km s _1 ) 
is the superbubble R2 and N18 (with the minimum veloc- 
ity of -328 km s" 1 ) is the SNR Sandage 8. It would be 
important to have kinematic data of the ionized gas over 
the whole dimensions of this galaxy in order to obtain its 
rotation curve. 

On the other hand, in Figure 14, the 21-cm isoveloc- 
ities from Lake & Skillman (1989) are sketched together 
with the radial velocity distribution of the Hn regions, 
obtained in this work. Their velocity resolution is 6.18 km 
s _1 and the linear resolution is approximately 210 pc. In 
this figure we can see that in general, the optical distribu- 
tion does not correspond to any particular kinematical H I 
feature. However, with respect to the density distribution 
(not shown), it has been noted that the optical extent of 
IC 1613 is located at the peak of the Hi density distri- 
bution and at the ridge of an Hi hole (Lake & Skillman 
1989). In order to compare the kinematics of the ionized 
gas with the H I kinematics we also require more observa- 
tions covering the whole optical extent of this galaxy. 

7. Conclusions 

We have measured - by means of a Fabry-Perot interfer- 
ometer - the kinematical parameters of the H n regions 
that we have detected within the 10' field of view of our 
observations (which includes most of the optical emission 
of IC 1613). We found that the Hn regions of IC 1613 
present similar integrated properties as well as Ha line 
profiles with simple structure that are fitted by single 
Gaussian curves. We obtained also the Ha fluxes and lu- 
minosities of these Hn regions. The luminosity function 
of these H II regions can be fitted to a power law function 
of slope a = —1.54 ±0.1 in agreement with the value given 



by HLG. We have also studied the diameter distribution 
and the cumulative diameter distribution of the Hil re- 
gions within our field of view. We obtained a characteristic 
diameter Do = 65.41 ± 3 pc which is slightly larger than 
the previous determination by HLG, probably because our 
field is smaller. Our characteristic diameter value does not 
seem to follow the correlation found by Hodge (1983) be- 
tween Dq and the absolute magnitude Mb for a sample of 
irregular galaxies. 

We have also found that the ISM of the entire galaxy is 
organized in a complex network of superbubbles. We have 
catalogued the superbubbles detected in our field of view. 
We have also calculated the 2-D porosity parameter, Qid-, 
in the field of view of our observations and seen that the 
Q2D value is intermediate between the observed values for 
the SMC and M31. 

Most of the superbubbles that we have detected show 
complex radial velocity profiles (some others are so faint 
that we cannot measure line-splitting in their velocity pro- 
files) . These kind of profiles along the same nebula are in- 
dicative of internal gas motions and expansions. We have 
determined the expansion velocities and other physical pa- 
rameters of the superbubbles found in this work. In ad- 
dition, we have used the photometric data published in a 
previous work (Paper I) and in Hodge (1978) in order to 
know the stellar content of the catalogued superbubbles. 
We find that almost every superbubble has an interior as- 
sociation containing massive stars, suggesting a physical 
link between them. Furthermore, we find that the dynami- 
cal time scales of the superbubbles are always shorter than 
the ages of the interior stellar associations. The mean tem- 
poral delay between those ages (about 10 Myr) suggests 
that only evolved members (i.e., the more massive stars) of 
these associations are responsible for the formation of the 
superbubbles, possibly when they evolve to supergiants or 
when they explode as supernovae. 

We have found that all the superbubbles with mea- 
sured expansion velocities have expansion velocities larger 
than 25 km s _1 . We used energy considerations in order to 
explain the origin of the superbubbles in terms of stellar 
winds and supernova explosions. We find that the super- 
bubbles with extreme expansion velocities (larger than 36 
km s _1 ) require interior stars with winds much more pow- 
erful than the average of B0-O9.5 supergiant stars if they 
are wind blown shells following the "standard model" . On 
the other hand, one or two supernova explosions could ex- 
plain the mechanical energies involved. However, no other 
traces of SN explosions (extended X-rays, non-thermal ra- 
dio emission) are known for these superbubbles except an 
enhancement of their [S 11] emission. 

We have obtained the size distribution and the cumu- 
lative diameter function of superbubbles detected in our 
field of view and compared these size distributions with 
size distributions of superbubbles in other galaxies. Our 
results are statistically poor. However, the value of the 
characteristic diameter we obtain, D = 91.84 ± 3 pc, is 
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Fig. 13. IC 1613 Ha observed radial velocity distribution (see text for details). The colors are related to velocities given in 
km s _1 . The scale and the orientation are indicated to the left 

in agreement with the value given by Price et al. (1990) 
for the nebulae in this galaxy. 

On the other hand, we have determined the global Ha 
velocity field of IC 1613, inside our 10' field of view. The 
velocity field is contaminated by the internal motions of 
the superbubbles of this galaxy. Nevertheless, a gradient 
of velocities is suggested across the SE-NW direction. In 
addition, we compare the mean velocities obtained from 
the H II regions with the kinematics of the H I gas but we 
do not find a clear agreement between them. 
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Table 1. IC 1613 properties 



Name IC 1613 

Type a IB(s)m 

R.A (1950) a 01h02mll.8s 

Dec (1950) +01d51m00s 

B magnitude a 9.88 

Angular size (arcmin) a 16.2 x 14.5 
Heliocentric radial velocity (km s^ 1 ) b -237 

Distance (kpc) c 725 

Inclination (degrees) 13 30 

P.A (degrees) b 50 



a NED 

b Optical heliocentric radial velocity from Martin (1998) 

c Adopted distance in this work, based on Freedman (1988a, 1988b) 



Table 2. Parameters of the FP Interferometric Observations 



Instrument 






PUMA 




Detector 






TEKTRONIX CCD 




FP Scanning Steps 






48 




Finesse 






24 




Filters 






Ha 


[Sir] 


Central Lambda (A) 






6570 


6720 


Interference order 






330 


322 


Free spectral range (km s 






847 


931 


Sampling resolution (km f 


5 1 /channel) 


18.9 


19.4 


Calibration line 






6562.78 H 


6717.04 Ne 


Binning 






2x2 


4x4 


Final pixel size (arcsec) 






1.17 


2.34 


Table 3. HLG HII regions 


taken for calibration 




Sandage a HLG b F c 


S d 


Counts Counts 




ID ID 






per sr f 




S2 13 220 


1.20 


27727 


1.51 




S18 31 131 


1.57 


15584 


2.00 




S6 36 108 


1.84 


26873 


4.72 




S16 46f 63 


5.90 


12398 


10.42 




46g 39 


3.47 


8721 


7.33 




59 45 


0.66 


7729 


1.12 




67c 105 


5.06 


15884 


7.16 




Sll 70 26 


2.09 


3489 


2.93 





a Sandage identification (Sandage 1971). 

b HLG identification (Hodge et al. 1990). 

c HLG published Ha fluxes, in units of 10 -15 erg cm -2 s _1 . 

d Surface brightness derived from HLG fluxes, in units of 10 -5 erg cm~ 2 s _1 sr _1 . 
° Counts (sky subtracted) in our image integrated over the areas defined by HLG. 
f In units of 10 12 counts sr _1 . 
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Table 4. Measured parameters for Hn regions in IC 1613 



Source a 


Sandage b 


HLG C 


D d 


Vi c 


v 2 




<T2 


log F(Ha) s 


log L(Ha) h 


Nl 


S2 


13 


127.6 


-248 




20.0 




-12.41 


37.34 


N2 




(25), 30 


53.3 


-242 




18.8 




-13.50 


36.25 


N3 




26 


78.1 


-254 




31.4 




-13.28 


36.47 


N4 




27,(33) 


47.2 


-242 




26.5 




-13.72 


36.03 


N5 


S12 


29,(35) 


80.5 


-247 




21.2 




-12.83 


36.91 


N6 


S18 


31,32 


99.0 


-247 




21.2 




-12.89 


36.86 


N7 




34 


34.2 


-244 




21.2 




-13.53 


36.22 


N8 


S6 


36 


57.5 


-240 




12.0 




-12.89 


36.86 


N9 


S3** 


37 


253.6 


-229 




18.8 




-12.14 


37.61 


N10 




38 


64.6 


-239 




18.8 




-13.36 


36.39 


Nil 


S7 


39,40 


56.6 


-239 




10.4 




-12.87 


36.88 


N12 




42,59 


117.3 


-229 




21.2 




-13.08 


36.67 


N13 




(44) 


29.3 


(-282) 




(5.9) 




-14.74 


35.01 


N14 


S17 


45 


90.8 


-236 




18.8 




-12.79 


36.96 


N15 


S15,S16 


46 


138.3 


-244 




22.2 




-12.14 


37.61 


N16 




47,57 


34.7 


-236 




17.6 




-13.27 


36.48 


N17 




48,58,61 


60.1 


-235 




16.3 




-13.09 


36.66 


N18 


S8* 


49 


34.5 


-239 


-335 


34.3 


86.1 


-12.54 


37.21 


N19 




50,41 


14.0 


-242 




17.6 




-13.61 


36.14 


N20 


S5 


51,52,53 


93.5 


-238 




16.3 




-13.05 


36.70 


N21 


S14 


55 


108.1 


-238 




18.8 




-12.63 


37.12 


N22 




56 


44.5 


-236 




22.2 




-13.28 


36.47 


N23 




(60) 


29.1 


-242 




22.2 




-14.51 


35.24 


N24 




63 


26.9 


-240 




24.4 




-13.94 


35.81 


N25 




64 


76.3 


-237 




24.4 




-12.87 


36.88 


N26 




66 


184.5 


-238 




21.2 




-12.46 


37.29 


N27 


S10,S13 


67 


198.1 


-248 


-185 


22.2 


18.9 


-11.98 


37.77 


N28 




69 


42.2 


-228 




29.5 




-13.63 


36.12 


N29 


Sll 


70,73 


42.7 


-235 




22.2 




-13.53 


36.22 


N30 




72 


99.0 


-228 




32.4 




-13.06 


36.69 



a Identification used in this work. 

b Sandage identification (Sandage 1971). The regions marked in boldface are known superbubbles previously catalogued as 
Hn regions. The asterisks indicate the regions hosting the supernova remnant (a single asterisk) and the Wolf-Rayet star 
(two asterisks). 

c HLG identification (Hodge et al. 1990). The regions marked in boldface are known superbubbles previously catalogued as 

Hll regions. The regions quoted between parentheses are objects with uncertain identification. 
d Equivalent diameter in parsecs (see text for definition). 
c Peak velocity in units of km s _1 . Parentheses indicate uncertain values. 
f Velocity dispersion in units of km s _1 . Parentheses indicate uncertain values. 
s Logarithm of the Ha flux in units of erg cm~ 2 s _1 . 
h Logarithm of the Ha luminosity in units of erg s _1 . 
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Table 5. Stellar content of the H II regions in IC 1613 


Source Known Related 


Stellar 


comments 


Stellar 


clusters e 


stars a object b 


association c 




association d 




Nl 






H5 


c32, c41 


N2 






H10 




N3 


G10 




Hll 




N4 


Gil 








N5 AM4 O-B 


G7 


embedded 


H10 




N6 


G10 


center 


Hll 


cl3 


N7 






H13 




N8 






H12 




N9 AM6 WO 










N10 


Gil 


center 


H13 


cl4 


Nil 


G14 


embedded 


H14 


c40,cl2 


N12 


G13 








N13 










N14 


G19 








N15 AM5 O-B 


G16, G17 


embedded 


H13 


ell 


N16 


G14 


embedded 


H14 




N17 


G18 


embedded 


H14 




N18 AM8 SNR 


G14 


embedded 


H14 


cl6 


N19 


G12,G14 


embedded 


H14 




N20 


G15 


center 


H15 


cl8 


N21 


G17, G21 


embedded 


H13 


cl5 


N22 










N23 










N24 


G24 


embedded 


H17 




N25 


G23, G24 


embedded 






N26 


G19 








N27 


G21, G25, G28 


embedded 


H17 




N28 


G28 


embedded 






N29 










N30 


G28 




H18 





a Stars with known spectra: WR candidates from Armandroff & Massey (1985) 

b Spectral types from Azzopardi et al. (1988). 

c Associations determined in Paper I. 

d Associations determined by Hodge (1978). 

e Stellar clusters determined in Hodge (1978) and in Paper I. 
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Table 6. Catalogue of Superbubbles in the galaxy IC 1613 



Source 


IDi a 
(ID 2 ) b 


axis c 
" x " 


Vi 

km s _1 


v 2 

km s _1 


log F(Ha) d 


log L(Ha)° 


Comments 


Rl 


GSla 
(S13) 


39 x 39 


-216 


-274 


-12.23 


37.51 


High [SII] emission 


R2 


GSlb 
(S10) 


49 x 33 


-244 


-147 


-11.98 


37.77 


High [SII] emission 


R3 


(HIV) (HHG 2) 


73 x 59 


-238 


-179 


-12.23 


37.52 


Faint [SII] emission 


R4 


GS2 

(HHG 1) 


66 x 50 


-277 


-213 


-12.18 


37.57 


Ring of HII regions 


R5 


GS4 
(S15) 


33 x 30 


-228 


-153 


-12.10 


37.65 


High [SII] emission 


R6 


GS5 

(S16+HIII) (HHG 4) 


64 x 52 


-240 


-190 


-12.64 


37.11 


Faint [SII] emission 


R7 


GS6 

(HIV) (HHG 5) 


47 x 38 


-234 


(-174) 


-12.72 


37.03 




R8 


(HVII) (HHG 3) 


63 x 49 


-259 


-208 


-13.08 


36.67 




R9 


(H5+S5+S7+S8) 


90 x 90 


-237 


-317 


(-12.31) 


37.44 


Ring of HII regions 


RIO 


(S6(contained)) 


37 x 37 


-229 


(-182) 


-12.71 


37.04 


[SII] detected 
S8 SNR f 


Rll 




77 x 77 


-205 




-13.19 


36.56 


Very faint 


R12 


(Hll(contained)) 


108 x 61 


-210 


-138 


-12.69 


37.06 




R13 


GS3 

(S12+S14) 


37 x 24 


-217 


-296 


-12.62 


37.13 


Expansion better 
seen in [SII] 


R14 


(HXII) 


85 x 77 


-234 




-13.39 


36.36 


Very faint 


R15 


(HIX) 


45 x 35 


-254 


-199 


- 13.20 


36.55 


Very faint 


R16 


(Dl) 


73 x 64 


-210 




-12.89 


36.85 


Very faint 


R17 


(HI contained) 


108 x 59 


-215 




-13.43 


36.32 





a Meaburn et al. (1988) identification. 

b Price et al. (1990) or HHG for Hunter et al. (1993) identifications. 

c Dimensions of the superbubble semiaxis in arcsecond units (see details in section 4.2). 

d Flux logarithm in units of erg cm -2 s _1 . 

e Luminosity logarithm in units of erg s _1 . 

f S8 SNR projected at the boundary of this superbubble. 
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Table 7. Stellar Content of the Superbubbles in the galaxy IC 1613 



Source 


Paper I 


N B 


Age 


Known 


Blue 


r 






Hodge (1978) 




(Myr) 


stars 


SG 






Rl 


part of G28 


3 


20 


no 


no 


-2.4 +/- 


0.3 




H18 














R2 


G25, part of G28 


17 


8 


no 


one 


-1.6 +/- 


0.2 




H17 








(B0) 






R3 


G29 


8 


20 


no 


no 


-2.1 +/- 


0.1 




H19 














R4 


parts of G16, G17, G21, G24, G25 


37 


8-10 


no 


no 


-1.6 +/- 


0.2 




H13, H17 (parts of) 














R5 


part of G16 


7 


10 


one 


one 


-1.7 +/- 


0.2 




H13 (part of) 






(OBem(?)) 


(S13) 






R6 


G5, G7, Gil 


36 


5 


one 


two 


-1.3 +/- 


0.3 




tji n ui o / t- 

HID, H13 (part of) 






(OB) 


(cB1.5, AOla) 






R7 


G9, part of Gil 


16 


10 


no 


no 


-1.7 +/- 


0.2 




mi (part or ) 














R8 


parts of , G9, G10 


5 


_ 


no 


one 


_ 














(B2I) 






R9 


G12, G14, G15, G18, part of G6, G13 


106 


5-10 


no 


no 


-2.0 +/- 


0.2 




H14, H15 














RIO 


G8, part of G12 


7 


5 


no 


one 


-1.1 +/- 


0.2 




part of H14 














Rll 


- 


— 


- 


— 


- 


- 




R12 


- 


- 


- 




- 






R13 


part of G21 


9 


8 


mi 


no 


-1.6 +/- 


0.2 




part of H17 














R14 




21 


16 


no 


no 








H6, H5 














R15 




36 


31 


no 


no 








H4 














R16 




8 


24 


no 


no 








H8 














R17 




25 


11 


no 


no 








H9 















22 
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Table 8. Physical Parameters of the Superbubbles in the galaxy IC 1613 



Source 


R 


n e (rms) 


Vexp 


K.E a 


t b 




pc 


cm~ 3 


km s _1 


ergs 


Myr 


Rl 


70 


3.8 


29 


1.88 


1.2 


R2 


74 


4.9 


49 


8.07 


0.7 


R3 


119 


1.8 


30 


4.64 


1.9 


R4 


105 


2.2 


32 


4.49 


1.6 


R5 


55 


6.3 


38 


2.56 


0.7 


R6 


105 


1.4 


25 


1.70 


2.0 


R7 


77 


1.9 


(30) 


(1.34) 


(1.2) 


R8 


101 


0.8 


26 


1.01 


1.9 


R9 


162 


(1.1) 








RIO 


66 


2.4 


(24) 


(0.69) 


(1.3) 


Rll 


138 


0.4 








R12 


152 


0.7 


36 


5.75 


2.0 


R13 


55 


3.6 


32 


1.04 


0.8 


R14 


146 


0.3 








R15 


72 


1.3 


28 


0.64 


1.2 


R16 


123 


0.7 








R17 


150 


0.3 









a Mechanical energy in units of 10 ergs. See section 5.2 for details. 
b Dynamical age in units of Myrs. See section 5.2 for details. 

Table 9. Comparison between the predicted wind power and the detected blue stars for the superbubbles of IC 1613 



Source 


Lvi/ a 


N B B 


Rl 


36 


3 


R2 


253 


17 


R3 


55 


8 


R4 


64 


37 


R5 


84 


7 


R6 


19 


36 


R7 


24 


16 


R8 


11 


5 


R9 




106 


R10 


12 


7 


Rll 






R12 


60 




R13 


29 


9 


R14 




21 


R15 


12 


36 


R16 




8 


R17 




25 



a Predicted wind power in units of 10 36 ergs s 1 . See section 5.4 for details. 
b Observed blue stars. See section 5.4 for details. 
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